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ABSTRACT 

t \ We investigate the consequences of a model of the assembly and growth of massive 

black holes (MBHs) from primordial seeds, remnants of the first generation of stars in 
a hierarchical structure formation scenario. Our model traces the build-up of MBHs 
from an early epoch, and follows the merger history of dark matter halos and their 
■ associated holes via Monte Carlo realizations of the merger hierarchy from early times 

I to the present time in a ACDM cosmology. The sequence of minor and major mergers 

(■f") , experienced by galactic halos in their hierarchical growth affects the merger history of 

MBHs embedded in their nuclei. So, if the formation route for the assembly of SMBHs 
<^> dates back to the early universe, a large number of BH interactions is inevitable. Binary 

• black holes coalescence timescales can be long enough for a third BH to fall in and 

interact with the central binary. These BH triple interactions lead typically to the final 
^ c"| expulsion of one of the three bodies and to the recoil of the binary. Also, asymmetric 

O ^ emission of gravitational waves in the last stages of the black hole merging can give a 

recoil velocity to the centre of mass of the coalescing binary. This scenario leads to the 
prediction of a population of intermediate mass BHs (IMBHs) wandering in galaxy 
halos at the present epoch. We compute the luminosity distribution produced by these 
I IMBHs accreting from their circumstellar medium. We find that in a Milky Way-sized 

galaxy they are unable to account for sources with luminosities > 10 39 erg/s unless 
they carry a baryonic remnant from which they are able to accrete for a long time. We 
>v> , also find that, for typical spiral galaxies, the bright end of the point source distribution 

j_j ■ correlates with the mass of the galaxy, and the most luminous sources are expected 

to be found in the disk. 

Key words: cosmology: theory - black holes physics - galaxies: kinematics and 
dynamics - X-rays 
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1 INTRODUCTION 

Observations of the centres of nearby galaxies indicate that 
most host a supermassive black hole (Magorrian et al. 1998; 
Ferrarese & Merritt 2000; Gebhardt et al. 2000); supermas- 
sive black holes (SMBHs) seem to be unequivocally linked 
to galactic bulges, their masses scaling with the bulge lumi- 
nosity and stellar velocity dispersion. In cold dark matter 
(CDM) 'bottom-up' cosmogonies, mergers of galaxies are a 
central part of the galaxy formation process in the hierar- 
chical structure formation picture, with major mergers (i.e. 
those between comparable-mass systems) expected to result 
in the formation of elliptical galaxies (Barnes 1988). 

Volonteri, Haardt & Madau 2003 (VHM03) and Volon- 
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teri, Madau, & Haardt 2003 (VMH03) have assessed a 
model for the assembly of SMBHS at the centre of galax- 
ies that trace their hierarchical build-up, assuming that 
the first 'seed' black holes (BHs) had intermediate masses, 
m. ~ 150 Mq, and formed in (mini)halos collapsing at 
2 ~ 20 from high-cr density fluctuations. These pregalactic 
holes evolve in a hierarchical fashion, following the merger 
history of their host halos. During a merger event BHs ap- 
proach each other owing to dynamical friction, and form a 
binary system. Stellar dynamical processes drive the binary 
to harden and eventually coalesce. 

The lifetime of BH binaries can be long enough 
(Begelman, Blandford & Rees; Quinlan & Hernquist 1997; 
Milosavljevic & Merritt 2001; Yu 2002) that following an- 
other galactic merger a third BH can fall in and disturb the 
evolution of the central system (VHM03). The three BHs 
are likely to undergo a complicated resonance scattering in- 
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teraction, leading to the final expulsion of one of the three 
bodies ('gravitational slingshot') and to the recoil of the bi- 
nary. Any slingshot, in addition, modifies the binding energy 
of the binary, typically creating more tightly bound systems 
(VHM03). 

Another interesting gravitational interaction between 
black holes happens during the last stage of coalescence, 
when the leading physical process for the binary evolution 
becomes the emission of gravitational waves. If the system 
is not symmetric (e.g. BHs have unequal masses or spins) 
there would be a recoil due to the non-zero net linear mo- 
mentum carried away by gravitational waves in the coales- 
cence ('gravitational rocket'). The coalesced binary would 
then be displaced from the galactic centre, leaving straight- 
away the host halo or sinking back to the centre owing to 
dynamical friction. 

Historically, the possibility of these gravitational inter- 
actions between BHs has been envisaged even before the 
CDM-hierarchical picture has become popular. Slingshots 
have been proposed, as a contrast to the unified -twin beam- 
theory (e.g. Blandford & Rees 1974) to explain the nature 
of double radio sources (Saslaw, Valtonen & Aarseth 1974; 
Mikkola & Valtonen 1990; Valtonen & Heinfimaki 2000): 
ejected BHs can produce extended radio lobes by interacting 
with the intergalactic medium in the areas of the lobes. 

Slingshots and rockets were both involved when trying 
to explain the nature of galactic dark halos as composed of 
BHs (Hut & Rees 1992; Xu & Ostriker 1994). Both processes 
basically give BHs a recoil velocity, that can even exceed the 
escape velocity from the host halo and spread BHs outside 
galactic nuclei (VHM03). So, the outlined scheme predicts, 
along nuclear SMBHs hosted in galaxy bulges, a number of 
wandering BHs that are the result of rockets, slingshots and 
mergers with a dynamical friction timescale longer than the 
Hubble time, so at z — the BHs are still on their way to 
the galactic centre. Note that most of the wandering BHs 
masses are in an intermediate range between stellar mass 
BHs and supermassive ones. 

In this paper, we improve the modelling of the dynami- 
cal interactions between galactic halos and between BHs and 
their host halos. In particular, we include the lengthening of 
the dynamical friction timescale of infalling satellites due 
to tidal interactions with the potential of the primary halo, 
and follow the orbital decay of every single satellite, deter- 
mining its position within the primary at any given time. 
We also trace step-by-step the trajectories of BHs ejected 
from galaxy centres following an energetic BH interaction; 
this enables us to keep track of the position and velocity of 
all wandering IMBHs as a function of cosmic time. 

The population of wandering IMBHs in local galaxies 
that our model predicts will have some observational con- 
sequences. In particular, due to accretion from their sur- 
rounding interstellar medium (ISM), these BHs will appear 
as bright X-ray sources. We compute here the luminosity 
distribution expected from such sources, and relate it to cur- 
rent observations. In particular, we discuss whether and un- 
der what conditions these IMBHs can make up a substantial 
fraction of a population of ultra- luminous sources (ULXs) , 
characterised by luminosities Lx > 10 39 erg/s, which has 
been recently identified in large Chandra surveys (e.g. Col- 
bert et al. 2004). The properties of this population as a 
whole (see Fabbiano 1988 for a review) hint that we might 



be actually observing an heterogeneous population, where 
super-Eddington luminosities could for example be due to a 
photon-bubble instability (Begelman 2002) or could be ap- 
parent due to strong beaming (King et al. 2001). A contribu- 
tion to the high-end tail of the X-ray point source luminosity 
distribution could also be due to young millisecond pulsars 
(Perna & Stella 2004). However, accreting IMBHs constitute 
an obvious candidate for the ULXs, as it has been often sug- 
gested (e.g. Colbert & Mushotzky 1999; Miller & Hamilton 
2002; Cropper et al. 2004). Part of our study here is there- 
fore a detailed investigation of the present-day observational 
properties of the population of IMBHs formed in mini-halos 
at high redshifts and evolved during the process of galaxy 
formation. 

The paper is organised as follows: in §2, we summarise 
the main features of our cosmological scenario of BH growth; 
the dynamical evolution of SMBH binaries is discussed in §3. 
The dynamical evolution of satellites BHs is studied in §4, 
while §5 describes the accretion model. §6 is dedicated to 
the observational consequences of these IMBHs in present- 
day galaxies and §7 discusses our results in the context of 
previous work in the field, and finally, in §8, we summarise 
our work. 



2 ASSEMBLY AND MERGERS OF SMBHs 

We briefly summarise here the main features of our scenario 
for the hierarchical growth of SMBHs in a ACDM cosmology 
(see VHM03 for a thorough discussion). In our model, pre- 
galactic 'seed' holes form with intermediate masses (m. = 
150 Mq) in (mini)halos collapsing at z = 20 from rare 3.5- 
<t peaks of the primordial density field (VHM03; Madau & 
Rees 2001). The assumed 'bias' assures that almost all halos 
above 10 11 Mq actually host a BH at all epochs. 

The merger history of dark matter halos and associated 
BHs is followed by cosmological Monte Carlo realizations of 
the merger hierarchy from early times until the present in a 
ACDM cosmology. The dynamical evolution of the BH pop- 
ulation is followed in detail with a semi-analytical technique. 
In this paper, we focus on the properties of two sets of ha- 
los with M = 2 x 10 12 M Q and M = 10 13 Mq respectively, 
creating 30 merger trees of the former and 10 of the latter 
and averaging the results. For the smaller halo, we will also 
consider a case in which seed holes are more numerous and 
populate the 3-<r peaks instead. 

In order to follow the history of halos and massive black 
holes, we adopt a two-component model for galaxy halos. 
The dark matter is distributed according to a NFW profile 
(Navarro et al. 1997). During the merger of two halo+BH 
systems of comparable masses, dynamical friction against 
the dark matter background drags in the satellite hole to- 
wards the centre of the newly merged system, leading to the 
formation of a bound BH binary. When two halos of mass M 
and M s merge, the 'satellite' (less massive) progenitor (mass 
M a ) is assumed to sink to the centre of the more massive 
pre-existing system on the dynamical friction (against the 
dark matter background) timescale, which depends on the 
orbital parameters of the infalling satellite, which we take 
from van den Bosch et al. (1999). 

The fate of the infalling halo has been investigated re- 
cently (Taffoni et al. 2003) in a NFW profile. Dynamical 
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ing on each BH individually. But as the binary separation 
decays, the effectiveness of dynamical friction slowly de- 
clines because distant encounters perturb only the binary 
centre of mass but not its semi-major axis. The BH pair 
then hardens via three-body interactions, i.e., by capturing 
and ejecting at much higher velocities the stars passing by 
within a distance of order the semi-major axis of the binary. 
The merger timescale is computed adopting a simple semi- 
analytical scheme that qualitatively reproduces the evolu- 
tion observed in N-body simulations (Merritt 2000; VHM03, 
VMH03). The binary evolution is slowed down due to the 
declining stellar density, with a hardening time that becomes 
increasingly long as the binary shrinks. If the hardening con- 
tinues sufficiently far, gravitational radiation losses finally 
take over, and the two BHs coalesce in less than a Hubble 
time. 

We have adopted a simplified description of the BH co- 
alescence process, ignoring for instance the depopulation of 
the loss cone (e.g. Yu 2002) and the loss of orbital angular 
momentum to a gaseous disc (Armitage & Natarajan 2002, 
Escala et al. 2004). 



Figure 1. Dynamical friction timescale of the halos hosting mas- 
sive black holes as a function of redshift ( left panel) and mass ratio 
(right panel). Note at low redshift (z < 5) the sharp division in 
timescale for major and minor mergers, the former characterised 
by timescales shorter than 2 Gyrs, the latter can be longer than 
the Hubble time. 

friction appears to be very efficient for 'major' mergers with 
(total) mass ratio of the progenitors, P — M s /M > 0.1. In 
each major merger the more massive hole accretes at the 
Eddington rate a gas mass that scales with the fifth power 
of the circular velocity of the host halo: the normalisation 
is fixed a posteriori in order to reproduce the observed local 
ibbh — cr» relation (Ferrarese 2002) and the observed lumi- 
nosity function of optically-selected quasars in the redshift 
range 1 < z < 5. 

Satellites of intermediate mass (O.OlMh < M Sj o < 
O.lMh) suffer severe mass losses by the tidal perturbations 
induced by the gravitational field of the primary halo. Tidal 
perturbations cause the progressive stripping of the satel- 
lite and, as the magnitude of the frictional drag depends on 
the mass of the satellite, this progressive mass loss increases 
the decay time. After a Hubble time satellites have typical 
masses ~ 1— 10% of their initial mass. The lightest satellites 
(P < 0.01) seem to be almost unaffected by orbital decay, so 
they survive and keep orbiting on rather circular, peripheral 
orbits. 

We adopt the relations suggested by Taffoni et al. for the 
orbital decay of satellites in a NFW halo. FigureQshows the 
evolution of the dynamical friction timescales with redshift 
and the dependence on the mass ratio for halos both hosting 
massive black holes. 



3 BLACK HOLE BINARY DYNAMICAL 
EVOLUTION 

The subsequent evolution of a BH binary is determined by 
the initial central stellar distribution. The binary will ini- 
tially shrink by dynamical friction from distant stars act- 



3.1 Gravitational Rocket 

When the binary semi-major axis is such that the timescale 
for hardening by stellar scatterings equals the timescale for 
shrinking by emission of gravitational waves, the latter pro- 
cess takes over. At this point, if the holes have unequal 
masses, or spin, the binary recoils due to the non-zero net 
linear momentum carried away by gravitational waves in the 
coalescence ('gravitational rocket'). 

The recoil velocity has still large uncertainties, and we 
adopt here the results from recent calculations by Favata 
et al. (2004). They estimate the recoil velocity as a func- 
tion of the binary mass ratio q = mz/mi (with BH masses 
mi 7112) and the spin parameter during final plunge and 
coalescence. We adopt here the expression they suggest for 
estimating 'upper limits' for the recoil, converting it into 
estimates of the binaries recoil velocity, i>cm, following the 
procedure described in Merritt et al. (2004). 

If the recoil velocity of the coalescing binary exceeds the 
escape speed v CS c, the holes will leave the galaxy altogether. 
If instead dcm < u csc , we assume that the binary is kicked 
out to a radius r max such that the total energy is conserved, 
i.e. we solve for r max the equation (mi +7TI2) Vq M = 2(\(f>o\ — 
|<?K r max)|), where <f>o is the galactic potential at the initial 
location of the binary and </>(r max ) is the potential at the 
apocentre. 



3.2 Gravitational slingshot 

The dynamical evolution of SMBH binaries may be dis- 
turbed by a third incoming BH, if another major merger 
takes place before the pre-existing binary has had time to 
coalesce or a wandering BH can reach the centre owing to 
dynamical friction (e.g. Hut & Rees 1992; Xu & Ostriker 
1994). If the incoming hole reaches the sphere of influence 
of the central binary, the three BHs are likely to undergo 
a complicated resonance scattering interaction, leading to 
the final expulsion of one, typically the lightest, of the three 
bodies (gravitational slingshot). 
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Typically an encounter with an intruder of mass mint 
smaller than both binary members leads to a scattering 
event, where the binary recoils by momentum conservation 
and the incoming lighter BH is ejected from the galaxy nu- 
cleus. By contrast, when the intruder is more massive than 
one or both binary components, the probability of an ex- 
change is extremely high: the incoming hole becomes the 
member of a new binary, and the lightest BH of the original 
pair gets ejected (Hills & Fullerton 1980). 

All triple interactions are followed along the merger 
tree, adopting the same scheme described in VHM03. 

In all cases we have estimated the recoil velocities 
conserving the binary energy and momentum and esti- 
mated the apocentre r max solving the equation (mi + 
m 2) v CM = 2(|^o| — |$(fmax)|)- Dynamical friction timescales 
have been calculated by the Chandrasekhar dynamical fric- 
tion timescale for circular orbits, so they are presumably 
upper limits to these timescales (thus giving an upper limit 
on the number of wandering BHs) , as the ejection probably 
occurs on eccentric orbits. 

In figure [5] the distribution of BH recoil velocities is 
plotted in units of the stellar velocity dispersion, ranging 
from 5kms -1 to 200 kms -1 depending on redshift and halo 
mass. As a rule of thumb, if the recoil velocity is smaller than 
the velocity dispersion, the BH remains in the proximity of 
the centre and the orbital decay is fast (Madau & Quataert 
2004). Figure [5] shows that the ejection of all three BHs is 
therefore not a likely outcome under the assumptions of our 
BH growth model (see VHM03 for a thorough discussion). 

Clearly, most of the single holes either escape their host 
halos or are slung to the periphery of the galaxy with con- 
sequently long dynamical friction timescales. Most of the 
binaries recoil instead within the core and fall back to the 
centre soon afterwards, with tdi < 0.01 Gyr. 

Note that triple interactions are rare events, in the 
whole history of a MW-sized halo with its satellites, on av- 
erage only ~ 3 slingshots happen, in the case of seed holes 
in 3.5-a peaks. 



4 DYNAMICAL EVOLUTION OF 
OFF-NUCLEAR BLACK HOLES 

The assumptions underlying our scenario lead to the pre- 
diction of a population of massive BHs wandering in galaxy 
halos, due to the processes described above. 

During the evolution and mass growth of DM halos, we 
track the position of all wandering BHs, following the mod- 
ifications of the galactic potential due to mergers. In case of 
minor mergers, the potential varies slightly, so the wandering 
BHs orbits evolve smoothly, while in case of major mergers 
the potential of both interacting systems change abruptly, 
causing larger modifications to the BH orbits. Although nu- 
merical simulations with resolution high enough to follow 
the dynamical evolution of intermediate mass BHs in galac- 
tic mergers are not available, we can rely on simulations 
studying the mixing of particles in major mergers (White 
1980, Nipoti et al 2003), assuming that during a merger a 
BH follows the evolution of other particles. On average par- 
ticles which were at a given radius (normalised at the half 
mass radius) in one of the galaxies before the merger, are 
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Figure 2. Distribution of BH recoil velocities (in units of the 
stellar velocity dispersion). Binaries, either following a triple in- 
teraction (upper left panel) or a rocket (upper right panel) , have 
a broad velocity distribution, while almost all single BHs have 
high recoil velocities, typically larger than the escape velocity. 



found at a similar radius (normalised to the new half mass 
radius) in the merger remnant. 

To track the orbital evolution of wandering BHs during 
galactic mergers we adopt the following scheme: 

• MINOR MERGERS: we let the wandering BHs in the 
more massive progenitor unperturbed, while the wandering 
BHs in the satellite galaxy are assumed to decay with their 
host. 

• MAJOR MERGERS: each wandering BH outside the 
core is placed in the new halo at a distance conserving the 
same ratio with the virial radius as it had in its parent sys- 
tem. If a wandering BH is instead inside the core, it is as- 
sumed to reach the centre of the newly formed halo within 
the same merging timescale of the halos (see section 2). 



5 DETECTING INTERMEDIATE-MASS 
BLACK HOLES IN LOCAL GALAXIES 

5.1 Naked wandering BHs 

As our first model of IMBHs in present-day galaxies, we 
assume that all BHs are naked, i.e. that the BHs left 
over by minor mergers have been completely stripped of 
their envelope. We embed all wandering BHs in a galactic 
disc/bulge/halo model (Springel & White 1999), and deter- 
mine the corresponding density of the medium in the sur- 
rounding of the BHs at their respective positions. We as- 
sume that a fraction ma of the initial halo mass is into a 
thin stellar disc with an exponential surface density, 



^(i^) = gexp(-|-)sech 2 (± 



(1) 



Dynamical evolution and observable signatures of IMBHs 5 



Here = m d M is the total mass of the disc and Rd is 
its scale radius. We adopt md = 0.05, while the disc scale 
length is determined through the fitting formula provided 
in Shen et al. (2003) for the Petrosian half-light (or mass) 
radius: 



R d = 0.17M^ 14 (1 + M d /3.98 x 10 10 ) 



(2) 



Most spiral galaxies seem to be consistent with a constant 
vertical scale length with a value of zo ^ 0.2i?d- Bulges are 
modelled with a spherical Hernquist (1990) profile of the 
form 



Pb(r) = 



2tt r(r b + r) 3 ' 



(3) 



In analogy to the treatment of the disc, we assume that the 
bulge mass is a fraction — 0.02 of the halo mass and that 
the bulge scale radius rb is a fraction /b =0.1 of that of the 
disc, i.e. rb = fbRd- We have then assumed that the gas 
traces the baryons, the gas fraction (of the total baryonic 
mass) is 20% in the disk and 4% in bulge and halo (Bell & 
de Jong 2000, Fukugita, Hogan & Peebles 1998). 

To model an elliptical galaxy (in the case M = 
10 13 M ) we have considered the Hernquist profile only, as- 
suming that the baryonic content of the galaxy is a fraction 
0.1 of the halo mass. For the elliptical, we have then assumed 
that the gas fraction is l%(Mathews & Brighenti 2003). 

5.2 Accretion model 

In this section, we discuss how we estimate the luminosities 
of off-nuclear BHs resulting from the process of accretion 
from the surrounding medium. 

Let p be the density of the medium in the surrounding 
of a BH, and let c s be its sound speed. The accretion rate 
onto a BH of mass Mbh can be estimated using the Bondi- 
Hoyle formula (Bondi & Hoyle 1944), 



M] 



Bondi — 



A 4tt G 2 Mbh P 

(V 2 +d) 3 /2 



(4) 



where v is the velocity of the BH with respect to the medium, 
and A is a parameter on the order of 1. The BH radial veloc- 
ities are computed by integrating the Chandrasekar formula 
(Binney & Tremaine 1987) as a function of the initial con- 
ditions, while the sound speed velocity, c s ~ 10(T/10 4 K) 1/ ' 2 
km/s, depends on the local temperature at the BH location, 
which is discussed later in this section. 

The bolometric luminosity of the BH can be written as 



Lboi = r)Mc 



(5) 



where r\ represents the fraction of the accreted mass that 
is radiated away. The nature of the accretion process, and 
the consequent value of r) in BHs is rather uncertain. Ac- 
tive Galactic nuclei (AGNs) accrete through accretion discs 
with an high efficiency (r) ~ 0.1). Their luminosities are 
close to their Eddington values. At the other end of the lu- 
minosity function, are the SMBHs at the centres of our own 
and nearby galaxies, whose luminosities can be as low as 
~ 10~ 9 — 10~ 8 of their Eddington values (e.g. Loewenstein 
et al. 2001). Low luminosities could be due either to low ef- 
ficiencies with M ~ M Bon di, or to relatively high efficiencies 
with M <C M Bon di (or also to a combination of both low ac- 
cretion rates and low efficiencies). Both scenarios have been 



suggested in the literature. Narayan & Yi (1995a) found 
that, when the accretion rate decreases below some criti- 
cal value rh ~ 0.1a 2 (a being the viscosity parameter in 
the disc), the thin disc switches to an Advection Dominated 
Accretion Flow, characterised by a low radiative efficiency, 
due to the inability of the electrons in the plasma to cool 
(see also Rees et al. 1982; Begelman & Celotti 2004). On 
the other hand, an accretion rate at values below the Bondi 
rate could be due to the presence of winds (or outflows; 
Narayan & Yi 1995b; Blandford & Begelman 1999). This 
scenario appears to be supported by observations of both 
accreting black holes (Loewenstein et al. 2001), as well as 
old, accreting neutron stars (Perna et al. 2003). 

In this work, when considering a model for the X-ray 
luminosity from accretion onto the IMBHs similar to that 
in SMBHs, we will adopt an empirical approach, and simply 
assume that the X-ray luminosity is a fraction r/x of the ac- 
cretion luminosity, with r/x calibrated on X-ray observations 
of SMBHs. This approach is similar to the one adopted by 
Agol & Kamionkowski (2002), in their calculations of the X- 
ray luminosity from isolated, stellar-mass BHs in the Milky 
Way (they parameterised their results in terms of the X-ray 
accretion efficiency). 

As a second model, we consider an optically thick, ge- 
ometrically thin accretion disc, as in the standard Shakura- 
Sunyaev (1973, SS) formulation. Radiation is emitted locally 
with a blackbody spectrum of temperature 



T(r) 



( 3GM BH M \ 
I 8nr 3 a ) 



1/4 



(t) 



1/2 



1/4 



(6) 



The emitted spectrum (per unit frequency) from the whole 
surface of the disc is given by 



2^ 



r dr 



e hv/kT{r) _ 1 



(7) 



We take Ri n = 3Rs and _R m ax = 1000i?s, where r s = 
2GMbh/c 2 is the Schwarzschild radius of the BH; this yields 
an efficiency rj ~ 0.06. As described in the previous sections, 
our simulations allow us to identify whether the BHs are in 
the halo, bulge or disc of the galaxy. For those in the halo 
and bulge, we assume that the temperature is the virial tem- 
perature of the accreting interstellar medium (ISM) . This is 
T vir = 5.8 x 10 5 K for the galaxy of mass 2 x 10 12 M Q , and 
T vir = 2.3 x 10 6 K for the galaxy of mass 10 13 M . In our 
Galaxy, we know that the disc is composed of several phases 
(e.g. Bland-Hawthorn & Reynolds 2000), from cold (T ~ 100 
K), to warm (T ~ 8000 K), to hot (T ~ 10 6 K). Modelling 
the various phases is beyond the scope of this paper essen- 
tially because, while some information on these phases is 
available for our Galaxy, no such studies have been possible 
for other galaxies. Moreover, due to the process of radia- 
tive feedback (Blaes, Warren & Madau 1995), the gas in the 
vicinity of the BH will be heated by the ionising radiation 
produced by accretion. The amount of heating and cooling 
depends on a number of factors, such as the density of the 
medium, its metallicity, the velocity of the BH, the accretion 
efficiency. For a medium made of H and He only, the ratio 
between the total heating rate to the initial energy density of 
the gas at the accretion radius is ~ 2 x 10 5 r;wJ 2 r 4 _1 (Agol & 
Kamionkowski 2002), where T 4 = T4IO 4 K, and v = u 4 o40 
km/s is the BH velocity. For BHs with velocities smaller 
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Figure 3. Mass function of wandering BHs at different redshifts, 
averaged over 30 Monte Carlo realizations of a M = 2 X 10 12 Mq 
halo assuming that seed BHs form in 3-cr density peaks at z=20. 
Upper panel: BHs in satellite remnants, lower panel: BHs dis- 
placed from the centre due to slingshots or rockets. Solid line: 
z=0, dashed line: z=2, dotted line: z=5. 



Figure 4. Mass function of wandering BHs at z=0, averaged over 
30 Monte Carlo realizations of a M = 2 X 10 12 M Q halo. Upper 
panel: seed BHs form in 3.5-cr peaks, lower panel: seed BHs form 
in 3-cr peaks. Solid line: total mass function, long dashed line: 
BHs displaced from the centre due to slingshots or rockets, short 
dashed line: BHs in satellite remnants. 



than a few tens of km/s, heating is important even for very 
low accretion efficiencies. Once H and He are ionized, the 
opacity drops substantially, and the heating rate becomes 
much longer than the accretion rate, leaving the gas at a 
temperature of around 2 x 10 4 K (Agol & Kamionkowski 
2002; Blaes, Warren & Madau 1995). We will take this to be 
the temperature of the accreting gas in the disc, assuming 
that, if there is a hot phase, this does not generally domi- 
nate. Note that, in those cases where the velocity of the BH 
is sufficiently high that dynamical heating is not important 
(and therefore the gas is actually cooler than assumed), the 
velocity of the BH dominates over that of the sound speed 
in the expression for the accretion rate, hence making the 
precise value of the temperature of the gas unimportant for 
our calculations. 

For any BH, the density of the accreting material is 
determined at its location in the galaxy according to the 
model described in §6.1. For satellites, we also consider the 
possibility that they might accrete from a remnant core as 
described in §5.4. 



6 RESULTS AND IMPLICATIONS FOR ULXS 

At a given time, a galaxy not only contains off-centre BHs 
slingshot or rocketed, but can also have BHs, possibly em- 
bedded in the remnant of their original host, still on their 
way to the centre following a minor merger. 

Figure |3| shows the evolution with redshift of the mass 
function of wandering BHs in a M = 2 x 10 12 M halo, 
assuming that seed holes form in 3-cr density peaks at z=20. 
Figure^Jshows instead the mass function of wandering BHs 
in a M = 2 x 10 12 Mq halo at z = 0, for the case of seed 



holes in 3.5-cr peaks and in 3-cr peaks. Note that most of the 
wandering BH masses are in an intermediate range between 
stellar mass BHs and supermassive ones. The upper end 
of the mass distribution of IMBHs, as well as their total 
number, is an increasing function of the galaxy mass. This 
can be seen by comparing Figure 4 (top panel) and Figure 
|S] which shows the mass function of wandering BHs in a 
M = 10 13 Mq halo at z = 0, for the case of seed holes in 
3.5-cr peaks. 

Figures [(J and |7| show the corresponding radial distribu- 
tions of these BHs. Typically, BHs in satellite remnants dwell 
in the outskirts of galaxies, while those caused by slingshots 
or rockets have a broader distribution and can be found even 
in the inner regions. Interestingly, there is not a correlation 
between BH masses and radial distance (Figure |HJ. Naively, 
more massive BHs would be expected in inner regions, due 
to the inverse dependence of orbital decay timescale on the 
mass. On the other hand, the shortness of the orbital decay 
timescale for the more massive BHs implies that very few of 
these are expected in the inner regions, because once they 
sink to the inner regions, they quickly sink further to the 
very center and merge with the central black hole or are re- 
ejected. Larger holes, however, appear only at later times, 
so that the time elapsed from their initial interaction with 
the galaxy is short. Figure |H| shows also the density of gas 
surrounding the same BHs. 

Resulting luminosities are shown in Figures 9, 10 and 
11. For the bolometric luminosities, we show results for an 
efficiency 77 = 0.1 (more typical of AGNs) and for rj = 10 -5 
(more representative of the dormant SMBHs in local galax- 
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Figure 5. Mass function of wandering BHs at z=0, averaged 
over 10 Monte Carlo realizations of a M = 10 13 Mehalo. Solid 
line: total mass function, long dashed line: BHs displaced from 
the centre due to slingshots or rockets, short dashed line: BHs in 
satellite remnants. 



Figure 7. Radial distribution of wandering BHs at z=0, averaged 
over 10 Monte Carlo realizations of a M = 10 13 Mq halo. Long 
dashed line: BHs displaced from the centre due to slingshots or 
rockets, short dashed line: BHs in satellite remnants. 
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Figure 6. Radial distribution of wandering BHs at z=0, averaged 
over 30 Monte Carlo realizations of a M = 2 X 10 12 Mq halo. 
Upper panel: seed BHs form in 3.5-cr peaks, lower panel: seed 
BHs form in 3-<r peaks. Long dashed line: BHs displaced from 
the centre due to slingshots or rockets, short dashed line: BHs in 
satellite remnants. 



ies), while X-ray luminosities are shown for the disc model 
and for an X-ray efficiency 1 7/x = 1CP 6 . 

For our galaxy model with mass M = 2 x 1O 12 M , 
and accretion from the local ISM, we find (see Figure 9) 
that there are no sources with bolometric luminosity larger 
than ~ 10 39 erg/s, even for a large efficiency of 0.1. Clearly, 
these sources cannot contribute to the ultra-luminous X-ray 
sources, with Lx > 10 39 erg/s, observed in local galaxies. 
We find the X-ray luminosity to be at most a few xlO 38 
erg/s (see Figure 10), if accretion proceeds at high efficiency 
through an SS disc. We need to stress that, although we 
find some very massive slung black holes at typical ISM 
densities (see Fig. 8), they have rather large velocities, and 
this suppresses the accretion rate and prevents these sources 
from exhibiting very large X-ray luminosities It should also 
be noted that, although the number of BHs is largest in 
the halo and the smallest in the disc, the bright end of the 
luminosity function is dominated by the BHs in the disc 
when accretion is solely from the ISM. This is due to the 
typically larger densities and lower temperatures of the disc. 
If the accretion mode of these BHs were similar to the one 
in the Galactic Centre, the bright end of X-ray luminosity 
would barely reach the ~ 10 34 erg/s level, several orders of 
magnitude below the Eddington limit for a solar mass BH. 

In our model in which local IMBHs are formed through 
mergers of smaller BHs during the galaxy assembling pro- 
cess, bigger galaxies are more likely to harbour more massive 
wandering BHs (compare Figures 3 and 4). Given the scaling 



1 The SMBH in the Galactic Centre has an X-ray efficiency of 
~ 4 X 10~ 8 - 2 X 10~ 6 in the 2-10 keV band (Baganoff et al. 
2001). For BHs at the other end of the mass spectrum (i.e. the 
stellar-mass BHs), Agol & Kamionkowski (2002) used a best guess 
estimate of 7)x ~ 10~ 5 . 
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Carlo realizations of a M = 2 X 10 12 Mq halo, with seed BHs in 
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Figure 10. Cumulative distribution function for the X-ray lu- 
minosity (0.1-10 keV band) of IMBHs accreting from the local 
ISM through a SS-disc (thick lines) and for an X-ray efficiency of 
10~ 6 (thin lines). The galaxy, of mass 2 X 10 12 Mq, is assumed to 
be a spiral, and the contribution from the disc, bulge and halo is 
separately shown. 




s 1 1 1 1 1 1 : 

- 3(7 satellite 




Disk 






Bulge : 






Halo 


n 








■i i < < 


M = 2xl0 12 M s : 
.1 1 1 1 



10» 10» 10 ; '« 10 : " 10™ 10 



1 34 1 35 1 38 1 37 1 3B 10 
L bol ( er § S "') 
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r\ = 0.1, while thin lines are computed with r\ = 10~ 5 . The contri- 
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is separately shown. The standard deviation among the various 
Monte Carlo realizations at the upper end tail of the distribution 
is 0.3 for the 3cr sling case, 0.22 for the 3cr satellites, and 0.18 for 
the 3.5cr sling. 



of the accretion rate with Mg H , this implies that, the larger 
the galaxy mass, the more luminous will be the sources as- 
sociated with accreting IMBHs. This is shown in Figure 11, 
where the bolometric luminosity function is computed for 
a galaxy of M — 1O 13 M0. Both a spiral and an elliptical 
galaxy models are considered. If the accretion efficiency is 
high (r) ~ 0.1), a spiral galaxy has a 10% probability of host- 
ing a source with L x > 10 40 erg/s. This is however not the 
case for an elliptical galaxy, where the lower ISM densities 
suppress the accretion emission. No ULXs due to IMBHs arc 
expected even in a very massive elliptical galaxy and for a 
high accretion efficiency. 



6.1 "Dressed" wandering BHs 

All the previous results have been produced under the as- 
sumption that the BHs left over by minor mergers have been 
completely stripped of their envelope and accrete from their 
surrounding ISM (naked BHs scenario). However, it is possi- 
ble that, during the merger process, BHs are not completely 
stripped of their surrounding core of material. Following 
tidal stripping, a bound core can remain, enveloping the 
central BH. Mayer & Wadsley (2003) have performed high- 
resolution N-Body/SPH simulations that follow the evolu- 
tion of small disc satellites falling into the potential of a 
much larger halo. Tidal shocks remove most of the outer 
dark and baryonic mass. In addition, ram pressure affects 
the gaseous component, stripping it even more severely. 
After a Hubble time, about 30% of the stars have been 
stripped, and 80 — 90% of the gas disc has been removed, 
leaving a system resembling a gas-poor dwarf spheroidal. A 
reasonable value for the gas density in satellite remnants, is 
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Figure 11. The upper panels show the bolometric luminosity 
distribution for a galaxy of mass M = 1O 13 M in a spiral (left) 
and elliptical (right) model. Thick lines are computed assuming 
an efficiency of 0.1, while thin lines of 10~ 5 . For the same galaxy 
models, the bottom panels show the X-ray luminosity (0.1-10 keV 
band) assuming that accretion proceeds through a SS-disc (thick 
lines) and for an X-ray efficiency of 10~ 6 (thin lines). Line types 
refer to the same galaxy components as in Figs. 9 and 10. The 
standard deviation among the various Monte Carlo realizations 
at the upper end tail of the distribution is 0.31 for the Spiral and 
0.48 for the Elliptical. 



the typical gas density p ~ 10~ 3 cm~ 3 , found in dwarf galax- 
ies of the local group (Grebel, Gallagher & Harbeck 2003) . 
The density scatter is large, varying by about 3 orders of 
magnitude from one system to another. An optimistic up- 
per limit for the gas densities was adopted by Islam et al. 
(2003b, c). Their estimate assumes that all the baryons in 
the satellite remnants are under the form of gas, and that 
the baryon content equals the cosmological baryon fraction. 
In this case, for a standard ACDM cosmology the gas density 
in the remnants is of order p~8x 10~ cm" . 

Given the uncertainties involved in the estimates of the 
density of these possible remnant cores, we prefer to just give 
some indicative results within the context of this scenario. In 
particular, we find that about 10% of the galaxies with mass 
~ 2 x 1O 12 M0 will have a source with bolometric luminosity 
of Lbol > 1.5 x 10 39 (j7/0.1)(p/0.1cm" 3 ) erg/s for the 3cr 
peaks, and L bol > 2 x 10 41 (r//0.1)(p/0.1 cm" 3 ) for the 3.5a 
peaks. It should also be noted that, if the BHs do accrete 
from a remnant core, then the luminosity function of the 
BHs will be independent of their location, and therefore the 
most luminous ones are likely to be found in the halo, where 
their number is the largest. 



Taylor & Silk 2003a,b,c). Although the hierarchical frame- 
work is similar this paper and the Islam et al. papers rely 
on different assumptions about the abundance of primordial 
seeds and about the dynamical and mass growth history of 
MBHs. First, their preferred model assumed seed holes, of 
m, = 260 M Q , in 3-a peaks at z — 24.6, leading to a num- 
ber of seeds almost 3 orders of magnitude larger than that 
in our reference model (3.5-cr peaks at z = 20). Second, their 
merger efficiency is much larger, since they assume that any 
MBH coming within 1% of the virial radius of the host halo 
actually merges with the central BH. Finally, they ignore the 
mass growth of MBHs due to accretion of gas, in a luminous 
phase as quasars, thus requiring a large number of MBH 
mergers in order to reproduce the observed local moH — c* 
relation. In our model, instead, it is gas accretion, follow- 
ing galactic major mergers, which determines the mass of 
SMBHs today, with mergers playing a secondary role, thus 
requiring a much smaller number of seeds to reproduce the 
observed SMBHs masses. 

As a result, our predictions for the number of X-ray 
sources in a galaxy of given mass generally fall below the 
predictions of Islam et al. (2003c), when similar estimates 
for the density and temperature of the core remnant are 
adopted . 



8 SUMMARY AND DISCUSSION 

In this paper we have studied the dynamical evolution with 
cosmic time of a population of IMBHs wandering in galactic 
halos, in the context of popular hierarchical structure forma- 
tion theories. We follow the merger history of dark matter 
halos and associated SMBHs via cosmological Monte Carlo 
realizations of the merger hierarchy from the end of the dark 
ages to the present in a ACDM cosmology. Massive black 
holes form as the end-product of the first generation of stars 
and then evolve in a hierarchical fashion following the fate 
of their host halos. Along the merger history of dark mat- 
ter halos, the BHs embedded in their inner regions sink to 
the centre of coalescing larger and larger galactic structures 
owing to dynamical friction, accrete a fraction of the gas in 
the merger remnant to become supermassive, and form a bi- 
nary system. Stellar dynamical processes drive the binary to 
harden and eventually coalesce. The dynamical evolution of 
SMBH binaries may be disturbed by a third incoming BH, 
if another major merger takes place before the pre-existing 
binary has had time to coalesce. The three BHs are likely 
to undergo a complicated resonance scattering interaction, 
leading to the recoil of the binary and of the single BH, 
thus creating a population of wandering BHs. Other coa- 
lescing binary BHs become wandering due to the non-zero 
net linear momentum carried away by gravitational waves 
if a binary is asymmetric. BHs retained within the poten- 
tial well of the host halo are then bound to slowly travel 
back towards the galactic centre under the effect of dynam- 
ical friction. The population of wandering BHs comprises 
also another class of holes, those originally hosted by halos 



7 COMPARISON WITH PREVIOUS WORK 

Our results are much more pessimistic than previous investi- 
gations involving IMBHs remnants of the first stars (Islam, 



2 Note that we consider the possibility of BHs accreting from a 
baryonic remnant only for the satellites, whose mass distribution 
has an upper end cutoff which is smaller than that of the sling 
BHs (see Fig. 3 and 4). 
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which are small compared to the primary halo. In this case 
the wandering holes can, or not, be still surrounded by the 
remnants of their old host. 

Our simulations have allowed us to estimate the number 
of black holes, their mass distribution and position within 
the galaxy. The total number of wandering IMBHs correlates 
with the mass of the galaxy at z — 0, and the BHs have typ- 
ically larger masses. As a result, the number of point sources 
due to IMBHs is also expected to correlate with the mass of 
the galaxy, and so the luminosity of the brightest sources. 
To zeroth order, no correlation is expected in our scenario 
of IMBHs formation and the star formation rate (SFR) in 
the galaxy at z — 0, as the mass function and positions of 
the BHs within the galaxy are essentially determined by its 
past merger history. However, if a galaxy has a high SFR, 
then it is likely to also have a larger filling fraction of dense 
molecular clouds, which would significantly enhance the ac- 
cretion rate if a BH happened to wander inside them. In 
the scenario recently proposed by Krolik (2004), BHs can 
be formed from primordial abundance gas or from stellar 
mergers 3 , and these mechanisms would make them to pref- 
erentially populate the disk of galaxies, therefore making 
them brighter, especially if they move for most of the time 
inside dense regions. 

In our scenario of BH formation in mini-halos and 
growth through subsequent galaxy merging, we find that, 
if the wandering IMBHs do not carry any substantial bary- 
onic core with them, then a Milky Way sized galaxy is not 
likely to harbour sources with L > 10 39 erg/s. On the other 
hand, about 10% of spiral galaxies with mass ~ 1O 13 M0 
would likely have a source with luminosity > 10 40 erg/s. The 
brightest sources are expected to be found in the disk, due 
to its typical lower temperature and larger density than the 
bulge and the halo. No ULXs are expected in ellipticals, even 
of large mass (~ 1O 13 M0), due to their average lower density 
and the lack of a cooler, dense disk which boosts the emis- 
sion. At sub-Eddington luminosities, the luminosity function 
of the naked IMBH population in a typical, MW-size galaxy 
is well below the observed point source population in Chan- 
dra surveys (Colbert et al. 2004). This is not surprising, as 
this point source population in the sub-Eddington range is 
dominated by high mass X-ray binaries (Grimm et al. 2004). 

Our results imply that, in order to produce ULXs from 
IMBHs in MW-type galaxies (and ellipticals for that mat- 
ter), it is necessary that the wandering IMBHs are not com- 
pletely stripped of their baryonic core, and this remnant 
material feeds the accretion at a rate higher than what 
would be allowed from the ISM alone. Even allowing an 
high efficiency of ~ 0.1, a density > 0.01cm -3 is required 
in order for a few percent of MW-type galaxies to host 
a source with luminosity > 10 40 erg/s. However, the den- 
sity cannot be significantly higher than that. In our models, 
the high tail of the mass distribution of IMBHs extends to 
~ 10 6 — 10 s Mq, and the observational lack of sources with 
luminosities > 10 44 erg/s implies that these IMBHs must 
accrete at sub-Eddington rates. 



3 For studies of stellar merging see also Portegies Zwart & McMil- 
lan (2002) and Gurkan, Rasio & Freitag (2004); note that the 
formation of IMBHs by runaway stellar merging is possible only 
in young, compact clusters. 



In this paper, however, we have investigated the lumi- 
nosity that IMBHs can attain in normal, quiescent galax- 
ies. Starbursting galaxies, such as the Antennae (Zezas & 
Fabbiano 2002) or the Cartwheel galaxy (Gao et al 2003) 
show much larger population of ULXs compared to quies- 
cent galaxies of the same size. In particular, in the Cartwheel 
galaxy these sources are found in the outer star-forming ring, 
believed to be created by radially expanding density waves 
caused by a plunging merger. Note that all these sources are 
at a distance larger than ~ 10 kpc from the centre of the 
galaxy, so in our scheme for a quiescent spiral, they would 
be found in the low-density, gas-poor galactic halo. The den- 
sity waves expanding outward, though, provide a temporary 
high density environment that may supply fuel for an IMBH 
to accrete and shine at high luminosities, higher also than 
those predicted for IMBHs travelling in galactic discs. 
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